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Abstract 

Under the assumption that cold dark matter and dark energy interact with each other through a 
small coupling term, Q, we constrain the parameter space of the equation of state w of those dark 
energy fields whose variation of the field since last scattering do not exceed Planck's mass. We use 
three parameterizations of w and two different expressions for Q. Our work extends previous ones. 
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I. INTRODUCTION 



The equation of state parameter of dark energy -the mysterious agent responsible for 
the current accelerated expansion of the Universe- ranks among the biggest unknowns in 
cosmology. It is usually written as w = p^/p^ and lies in the range (—1,-1/3) if the dark 
energy is quintessence, and (—00, —1) if it is of phantom type; in either case it may run with 
expansion. By contrast, in the very particular (but observationally favored) instance of a 
cosmological constant w stays fixed at —1. It is then apparent that some knowledge on the 
nature of dark energy can be gained by setting limits on w. 

Observational constraints on this quantity from supernovae type la (SN la), tempera- 
ture anisotropies of the cosmic microwave background radiation (CMB), baryon acoustic 
oscillations (BAO), etc, suggest that nowadays w is not far from the cosmological constant 
value. However, all these studies rest on some or other world model and on inevitable priors 
imposed upon the parameters space to ease the data analysis. 

Recently, bounds on w have been derived from the very reasonable requirement that the 
variation experienced by the dark energy scalar field <fi (regardless it may be phantom or 
quintessence) from any redshift, z, within the classical expansion era, till now should not 
exceed Planck's mass (see Huang jl], and Saridakis [2]). This is to say 



\A<p(z)\/M P <l, (1) 



where -as can be readily deduced- 



ia*mi_ 1 r\ dt= r vwr^sm^. (2) 



M P M P J t J l + x 

Huang's analysis confines itself to quintessence fields with various w parameteriza- 
tions (viz: w constant, a linear function of redshift, and the Chevallier-Polarski-Linder 
parametrization Saridakis considers phantom fields only and uses identical pa- 

rameterizations than Huang plus a further one in which w depends linearly on the logarithm 
of the scale factor 2)]. 

The bound (0Q) looks a rather natural condition and persuasive arguments have been 



advanced in its favor 



5H7(]. Moreover, we have found a further, solid, motivation for it. 
Namely: if condition ([1]) is violated, then dark energy dominates at very early times which 
cannot be reconciled with observation (see Fig. H] and "forbidden" regions on the top right 
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corner of the panels of figures HH [H and [HI below). In particular, as noted by Bean et 



al. [8], the fractional density of dark energy cannot exceed 5% and 39% at the primeval 
nucleosynthesis and recombination epochs, respectively, -we shall elaborate on this in a 
future publication. Nevertheless, we refrain ourselves from imposing ([T]) on every dark 
energy field. Its seems safer to demand that \A<p\/Mp be not much larger than order unity. 
However, the latter requirement cannot be easily implemented and, at any rate, it would 
translate on rather loose constraints on w. Therefore, we choose to circumscribe our analysis 
to dark energy fields that fulfill condition (TjQ). 

In Huang's as well as in Saridakis' work matter and dark energy interact with each other 
only gravitationally. Certainly this is a reasonable assumption but nonetheless minimalist 
since, while dark energy interactions with baryons are severely restricted by solar system 



experiments [9|, there is nothing against a possible coupling (interaction) with dark matter. 
On the contrary, a transfer from dark energy to dark matter alleviates the cosmic coincidence 
problem (the observational fact that both energy densities are comparable today [l0|) and 
may solve it in full [Uj. Besides, as revealed by optical, x-ray and weak-lensing data, the 
internal dynamics of 33 relaxed galaxy clusters seems to favor such interaction [12| . Likewise, 
a recent analysis using the 397 SN la of the "Constitution" set [131 ] . BAO and CMB data 



adds significant weight to the existence o 
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the coupling [yj]. In fact, the subject has evolved 



16] and references therein; a recent review can be 



into a field of active research -see e.g. 
found in [l7j. 

In the presence of a coupling, say Q, the continuity equations for the three main contrib- 
utors to the present cosmic budget adopt the form, 

p b + 3Hp b = , 
p c + 3Hp c = Q , 

p + 3H{l + w)p <t> = -Q ) (3) 

where the subscripts, b, c, and <ft stands for baryons, cold dark matter, and dark energy, 
respectively. This paper aims to constrain different parameterizations of w under diverse 
expressions for Q. Regrettably, guidelines about the latter are rather loose: it must be 
positive-semidefinite and small. If Q were negative, the energy densities could become 
negative, and the second law of thermodynamics get violated 18] . If it were large, dark 



energy would have dominated the expansion at early times and possibly not today. Thus, 



for the time being, we must content ourselves with guessing plausible expressions for Q just 
on phenomenological bases. 

In view of the above equations the coupling must be a function of Hp c and Hp^. By 
power law expanding it up to first order on these quantities we write Q(Hp c , Hp^) ~ e c Hp c + 
e^Hp^, where both e coefficients ought to be non-negative and small (i.e., not larger than, 
say, 10 _1 ). To simplify the analysis we will consider in turn that one of the two coefficients 
vanishes. So, we will take 

Q = 3 eHptp , and Q = 3 eHp c , (4) 

-the factor 3 being introduced for mathematical convenience. 

At this point one may object that if Q is small, no significant difference with the results 
of Huang [1] and Saridakis [2] should be expected. However, as we shall see, this is not 
;he case. A priori we can argue in favor of some non small departure from the findings of 
l| and [2I by considering the ratio between Q and the second term in last equation of 
Indeed, dividing the right hand side of (jUl) by the second term of the said equation yields 
e/(l + w). In general, this ratio cannot be neglected and it can be large in absolute value 
when w is close to —1. Likewise, using instead (JU2) we have ep c /(l + w)p ( f > . Its absolute 
value can be of order unity, or even larger, for extended periods of the cosmic history. 

Although expressions (jlj) were proposed on phenomenological grounds they can be ob- 
tained in scalar-tensor gravity from the action (in the Einstein frame) [3] 

S = j d'x^-g j-l^R-I^d'V + -^L c (CdC, X - l 9a b ) + L b {i,d^ X ~ l 9a^ ■ 

(5) 

Here R denotes the Ricci scalar, L c and ( stand for the dark matter Lagrangian and the 
dark matter degrees of freedom, respectively (corresponding meanings have Lf, and £); on the 
other hand, x(0) couples dark matter with the dark energy field <fi. From ([5]) the interaction 
term between dark matter and dark energy can be expressed as 

d(mx(a)) 



Q = H Pc 



(6) 



din a 

where x( a ) = x( a )^ 3u,c1 ^ 2 with w c the equation of state parameter of dark matter (zero in 
our case). By choosing 



x (a) = xo exp 



din a 



Pc 



(7) 
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the first (second) expressions for Q in PJ follows after using (jSJ) and setting e c = (e^ = 0). 



20| and |21] 



Similar expressions for Q have been obtained from the above action in Refs. 

In this paper we shall assume a spatially flat Friedmann-Robertson- Walker universe with 
present fractional densities of baryons and cold dark mater VL^q = 0.04, and Q c q = 0.24, 
respectively, throughout. As usual, a zero subindex means the current value of the corre- 
sponding quantity; likewise we normalize the scale factor of the metric by setting ao = 1. 

II. INTERACTION TERM PROPORTIONAL TO THE DARK ENERGY DEN- 
SITY 

In this section we consider Q = 3 eHp^ alongside different expressions for the equation of 
state parameter. 

A. Constant w 

By plugging Q = 3-f/ep^, into (jHJ), assuming w = Wq, and integrating we get the energy 
densities dependence on the scale factor 



p&aT* [1 - a - 3 < tt0+e >] , 



Pb 


= Pboa~ 


-3 

i 


Pc 


= Pcoa' 


- 3 + 


Pa 


= Pso Q> 


-3(H 



w 

-w +t) ^ ^ 

We then introduce = p^,/{pb + Pc + p<f>) in (EJ) and numerically determine \A(f>\/Mp in the 
redshift interval between the last scattering (zi s = 1089) and now (z = 0) -bear in mind that 
1 + z = a' 1 . Figure [1] shows \A<f>\/Mp vs. wq for different values of the parameter e. When 
— 1 < w < — 1/3 (left panel), for every e value there is a range of Wq that violates condition 
(0Q); the bigger the interaction strength, the bigger the \wq\ that violates the condition. By 
contrast, when wq < —1 (right panel) the bound (CQ) is fulfilled no matter the value of e. 
This could have been anticipated since, in this simple case, Eq. (J3J2) can be written as 
+ 3H(1 + w e ff) pj, = 0, i.e., as though it were no interaction but with w replaced by the 
effective equation of state parameter to e // = w + e. At any rate, our results corroborate and 
extend Huang's and Saridakis'. 
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FIG. 1: \A(p\/Mp vs. wo for Q given by (jUl) an d constant w evaluated in the redshift range 
(0,1089). Left Panel: Quintessence solutions (—1 < < —1/3). Right panel: Phantom so- 
lutions (wq < —1). In both panels solid, dotted, dashed, and dot-dashed lines correspond to 
e = 0, 10 -4 , 10 -2 , and 0.1, respectively. (Solid and dotted lines practically overlap). Condition ([1]) 
for e = 0, 10" 4 ,10~ 2 , and 0.1 implies w < -0.738,-0.738,-0.745, and -0.810, respectively. In 
drawing these and all subsequent figures we assumed f2&o = 0.04, and f2 c o = 0.24. 

B. Chevallier-Polarski-Linder parametrization 

Except when dark energy is given by the quantum vacuum, there is no compelling mo- 
tivation to assume w constant for the whole cosmic evolution. However, its simplest gen- 
eralization in terms of redshift, w(z) = Wo + W\Z, is not compatible with observation for 
it diverges as z — > oo. This prompted the introduction of the more suitable expression 
w(z) = wq + W\ or, equivalently, in terms of the scale factor 

w(a) = w + u>i(l — a) (9) 



by Chevallier and Polarski 



3| (later popularized by Linder [4]) which does not suffer from 
that drawback and behaves nearly linear in z at low redshifts. (Clearly, w = w(z = 0) and 

Wi = dw(z)/dz\ z=0 ). 

Now the dark energy density integrates to = p^ a~ 3 ( 1+Wo+Wl+e ^ exp [3wi (a — 1)]. By 
contrast, the energy density of cold dark matter has no analytical expression but it can be 
found by numerical integration of Eq. ([3j2). For the sake of illustration we present in Fig. 
[2] the evolution of its fractional density for a certain choice of parameters. Obviously, the 
expression for the baryon density does not vary. 




FIG. 2: Evolution of the fractional density of cold dark matter O c = p c /(pb + Pc + p<j>), since last 
scattering till now, in terms of the scale factor for the Chevallier-Polarski-Linder parametrization 
([9]) with wo = — 1 and W\ = 0.5. Solid, dotted, dashed, and dot-dashed lines correspond to 
e = 0, 10 -4 , 10~ 2 , and 0.1, respectively. The two first lines are undistinguishable from one another. 

Proceeding as before, we numerically evaluate \A<p\/Mp in terms of the parameters e, Wo 
and W\. Figure [3] shows w vs. W\ for different choices of e. Points in the plane (wi,w ) 
satisfying simultaneously wo+wi > —1 and wo > —1 fulfill w(a) > —1 for all a > 1 -i.e., they 
correspond to quintessence models that never evolve into phantom. For small e condition 
(|T]) is violated in a section of the quintessence region, but only there. However, when e 
takes moderate values, the condition gets also violated by some models that evolved from 
phantom to quintessence (some few models in the mixed region I) and from some models 
that evolved in the opposite sense (some models in mixed region II), see right bottom panel. 
Models that always stay phantom respect condition ([1]) for all e at any redshift z > 0. 
The triangular region wq + w\ > at the top right corner of each panel is observationally 
forbidden since points lying there correspond to models that feature dark energy dominance 
at high redshifts. It is noteworthy that it entirely falls in the region that violates condition 
([1]). Likewise when e is about 0.029 or larger, an unphysical region of negative dark matter 
density develops -again in the region that violates (JT])-, see right bottom panel. It is also 
visible that A0 slowly augments with e. 
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e-o. e- 0.0001 




FIG. 3: In the light-gray sections of the four panels condition (pQ) is met when use of both (jUl) anci 
the Chevallier-Polarski-Linder parametrization is made. In the dark-gray sections, \A(p\/Mp < 
0.7; they shrink with increasing e. The regions marked "Ever Quintessence" correspond to models 
such that wo + w\ > — 1 and wq > — 1 for all z > 0. Similarly, the regions marked "Ever Phantom" 
correspond to models such that w(a) < —1 for all a < 1. In the mixed region I of each panel 
models transit from phantom to quintessence as the Universe expands. In the mixed region II they 
transit in the opposite sense. For wo > —0.46 the expansion is decelerated at a = 1, as indicated 
in each panel, see the text. 

Also marked in each panel is the maximum wo for which there is acceleration at z = 0. 
This value, wq = —0.46 , readily follows by setting the present deceleration parameter, 
q = -d/(aH 2 )\ z=0 = \ [1 + 3w (1 - Q b0 - Q c0 )}, to zero. 

As we have checked, for w values on the border of the forbidden region the ratio A<fi/Mp 
is, at high redshifts (say, z > 20), consistently larger than 2 and can be as large as 7 -see 
Fig. HJ The fact that w values in the forbidden region correspond to dark energy dominance 
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FIG. 4: Acf)/Mp in terms of the normalized scale factor. In both panels e = 0.01 and the w values 
correspond to points on the straight border line of the forbidden region of the left bottom panel of 
Fig. [3J Left panel: wo = —0.46. Right panel: wo = —1.2. 

at early times and this is incompatible with observation (e.g., the primeval nucleosynthesis 
scenario and cosmic background radiation spectrum would be very different from the one 
admitted today, and galaxies could not have formed), suggests that, in general, scalar fields 
respect the bound ([I]). This also holds true for the forbidden regions of Figs. |6l |9j and [PT1 
below. 



C. Barboza-Alcaniz parametrization 

As readily noted, Chevallier-Polarski-Linder's parametrization ([9]) implies that w(z) di- 
verges as z —7- — 1 (i.e., in the far future). To avoid this unpleasant feature Barboza and 
Alcaniz proposed the ansatz w(z) = w + w\ Vffi or, equivalently, 

w{a) = Wo + Wl l-la + 2a* > (10) 
which ensures that w(z) stays bounded in the whole interval — 1 < z < oo aside from 



behaving linearly in redshift for \z\ <C 1 



22] 



Assuming this novel parametrization alongside the interaction term (jUl), we get = 
P<po a~ 3 ( 1+Wo+Wl+€ \l — 2a + 2a 2 ) 3wi / 2 . Again, p c must be calculated numerically. Figure E] 
illustrates the behavior of Q c . 

Proceeding as before, we numerically assess \A<p\/Mp in the redshift interval (0,zi s ). 
Figure |6] depicts w vs. Wi for various choices of e. If w\ > 0, quintessence models lie in the 
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FIG. 5: Evolution of the fractional density of cold dark matter tt c = p c / (pb+ Pc+ P</>) m terms of the 
normalized scale factor using the parametrization of w(a) of Barboza and Alcaniz (|10p with wq = — 1 
and w\ = —1.5. Solid, dotted, dashed, and dot-dashed lines correspond to e = 0, 10~ 4 , 10~ 2 , and 
0.1, respectively; the three first practically overlap each other. 

region given by — 1 < wq — 0.21 w\ and wq + 1.21 w\ < 1; if w\ < 0, they lie in the region 
— 1 < Wq + 1.21 W\ and Wq — 0.21 w\ < 1. Unlike the preceding section, in the mixed regions 
II of the parameter space, phantom models (that started as quintessence at high redshift) 
violate condition ([I]) for whatever e. As before, the triangular regions given by Wq + Wi > 
are observationally discarded as they correspond to dark energy dominance at early times. 
These ones are much wider than in the previous case. Again a region of negative cold dark 
matter density develops, this time for e values in excess of 0.015 -see right bottom panel. As 
it is apparent, the interaction has a noticeably effect on the variation of <ft i n the sense that 
A(p moderately augments with e. Comparison of Figs. |3] and [6] reveals that the interaction 
induces a bigger evolution of when w(z) follows Barboza- Alcaniz's parametrization than 
when it follows the parametrization of Chevallier-Polarski-Linder. 

III. INTERACTION TERM PROPORTIONAL TO THE COLD DARK MATTER 
DENSITY 

In this section we take up Q = 3 eHp c for the coupling term and resort to the expressions 
for the equation of state parameter considered in the previous section. 
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FIG. 6: Same as in Fig. [3] except that w(a) obeys the parametrization of Barboza and Alcaniz, 
Eq. (TO]). 

A. Constant w 

For w = Wq we readily find 

Pc = p c0 ^ 3(M , 

P* = P*o «- 3(1+, " 0) + Pco -4- [ a ~ 3il+W0) ~ «" 3(1 " £) ] • ("J 

w + e 

Fig. [7] shows \A<p\/Mp vs. Wq for different values of e. For quintessence fields there is 
a range of il>o that violates condition (CQ) regardless the value of e. For phantom fields the 
said condition only holds when e is very small. This contrast with non-interacting phantom 
models which fulfill \A<j)\/Mp < 1 irrespective of w (see Fig.l of Ref. Q). 



II 



-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -2.5 -2 -1.5 -1 

FIG. 7: \A(p\/Mp vs. wo for Q given by dH2) and constant w evaluated in the redshift 
range (0,1089). Left Panel: Quintessence solutions (—1 < w$ < —1/3). Right panel: Phan- 
tom solutions (wo < — 1). In both panels solid, dotted, dashed, and dot-dashed lines cor- 
respond to e = 0, 10 _4 ,10 -2 , and 0.1, respectively. Condition ^Q) implies w$ < —0.738 for 
e = and w < -0.745 for e = 10~ 4 ; for e = 10~ 2 it implies -1.299 < w < -0.841, and 
-1.071 < w < -0.939 for e = 0.1. 
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FIG. 8: Evolution of the fractional density of dark energy for the Chevallier-Polarski-Linder 
parametrization, Eq. ([9]), with wo = —1 and w\ = 0.5. Solid, dotted, dashed, and dot-dashed 
lines correspond to e = 0, 10~ 4 , 10~ 2 , and 0.1, respectively. 

B. Chevallier-Polarski-Linder parametrization 

In this case p c = p c o a~ 3 ^~ e ^ while the expression for must be found numerically. 
Figure |8] instances the evolution of Vt^ in terms of the normalized scale factor. 
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FIG. 9: Same as Fig. [3] except that the interaction term is now given by Eq. (JH2) and that in the 
right bottom panel only the \A(j)\/Mp < 1 region is depicted. 

Dark energy models in Fig. |9]such that Wq+W\ > — 1 and Wq > — 1 present a quintessence 
behavior during its whole evolution. For small e only models in this region can violate 
condition (pQ). For larger e values, models in the other regions (bottom panels) can also 
violate it. As it is manifest, the area of the offending region augments with e. As in 
the preceding cases, models in the upper right triangle of each panel are observationally 
discarded. 



C. Barboza-Alcaniz parametrization 

In this instance p c bears the same expression in terms of the scale factor than in the 
previous case. Figure ITU1 illustrates the dependence of on a. 

After having numerically evaluated \A<p\/Mp, Fig. (TT] displays u^o vs. W\ for different 
choices of e in the redshift interval (0,^ s ). If w\ > 0, ever quintessence models belong to 
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FIG. 10: Evolution of the fractional density of dark energy for the Barboza-Alcaniz parametriza- 
tion, Eq. (fTUj) . with wo = — 1 and w\ = 0.5. Solid, dotted, dashed, and dot-dashed lines correspond 
to e = 0, 10~ 4 , 10~ 2 , and 0.1, respectively. 
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" 1 "l 
£ = 0.01 £ = 0.1 




FIG. 11: Same as in Fig. [9] except that w(a) obeys the parametrization of Barboza and Alcaniz, 
Eq. (HOD. 
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the region given by —1 < wq — 0.21 wi and wq + 1.21 w\ < 1; if w\ < 0, they lie in the 
region — 1 < wq + 1.21 W\ and wq — 0.21 w\ < 1. Even for small e, models in this region 
as well as models in the mixed region II can breach condition 0JJ) (top panels). For larger 
e values, models in the other regions (bottom panels) can also breach it. As it is apparent 
from the right bottom panel (e = 0.1), models that always stay phantom are restricted to a 
very narrow region close to the line Wq = — 1 — 1.21w\. 



IV. DISCUSSION 



Motivated by the reasonable assumption that the variation A(f> of the field driving the 
present phase of cosmic accelerated expansion should not exceed Planck's mass we have nu- 
merically calculated the said evolution since last scattering till now for different expressions 
of the equation of state parameter w(a) and two couplings, Q, between dark matter and 
dark energy. This constrains the parameter space (w±, Wo) as shown in the six cases studied. 
Figures (TJ El El El and [TTJ explicitly exhibit the regions of the parameter space (w%, w ) 
ihat fulfill the bound \A<p\/Mp < 1. Our analysis extends those of Huang [l] and Saridakis 
^ who assumed that dark matter and dark energy evolved separately and did not consider 



Barboza and Alcaniz's parametrization 



22]. 



The main results of this work can be summarized as follows: 
(i) When the interaction obeys (jUl), models that always stay phantom respect condition 0JJ) 
at any redshift and for any e value, (ii) Likewise, depending on the values taken by wo and 
Wi models that always stay as quintessence as well as models that evolve from quintessence 
to phantom and models that evolve in the opposite sense can violate the said condition. (Hi) 
When Q obeys OH 2) and w(a) is given by the Chevallier-Polarski-Linder parametrization 
(E]), for small e, all models that ever stay phantom and all models in the mixed regions 
I and II respect bound ([T]); -top panels of Figs. El However, for not so small e values, 
models in these regions may also violate ([T|) -bottom panels of Fig. El (iv) When Q is given 
by 0H2) and w(a) obeys the Barboza- Alcaniz parametrization (flOl) . for small e, all models 
that ever stay phantom and models in the mixed region I satisfy bound ([1]) -top panels of 
Fig. [HI By contrast, for not so small e values, models in these regions may also violate 
01]) -bottom panels of Fig. [TTJ (v) The dark energy field <fi appears to evolve faster when: 
(a) the coupling Q depends explicitly on p c than when it does on p^, and (6) when w(a) 
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obeys the parametrization of Barboza-Alcaniz than when it does that of Chevallier-Polarski- 
Linder. Point (a) is in line with Caldera-Cabral et al. result that interacting models with 
Q = SeHp^ "work better" than models with Q = 3eHp c Further, these models may 



help explain the non-vanishing temperature (about 0.6 Kelvin |23|) of sterile neutrinos 24]. 

Our finding that | A0| increases with e strengthens the view expressed in the Introduction 
that Q must be small. Otherwise (except for flUl) and w = constant < —1, right panel of 
Fig. [1]) condition ([1]) would be violated for whatever w not far from —1. 

We have confined ourselves to models in which the interaction term is proportional to the 
Hubble factor (Eqs. 01])) whence our results apply to them only. Other interaction terms can 



be found in the literature, among others, Q oc cf> p c [25], and Q oc (p c + p^) -see e.g. [16l. |26|. 



However, the former class of models present the drawback of being unable to simultaneous 
lead to a correct sequence of cosmic eras (radiation, matter, and dark energy) and solve the 
coincidence problem j^, while the models we have considered do not suffer from that 28]. 



As for the latter class of models, they generally lead to negative energy densities, either of 



dark matter or dark energy, at early or late times 26]. Still, for some specific values of the 
parameters entering the interaction these models are free of that problem. 

As is well known, due to quantum instabilities phantom fields may find no place in Nature 
-see e.g. 29) . However, some phantom models based on low-energy effective string theory 



may not suffer from such pathologies 



Moreover, 



3l| 



301 ] whence this issue is not settled as yet. 
higher derivative terms in the Lagrangian may render phantom models stable c 

Before closing, it is sobering to recall that variable dark energy fields are afflicted by the 
problem of their small mass (m^ ~ 10~ 33 eV). Due to this unpleasant feature they are looked 
upon as not more natural than the cosmo logical constant, despite the enormous fine tuning 
scale of the latter. Admittedly, this is an unsolved problem that affects most (if not all) 
candidates of evolving dark energy. However, recently, a supergravity based approach that 



might point to the solution has been proposed 32]. It provides a very small mass for the 



field and a seemingly natural mechanism for a weak coupling between the field and matter. 
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